Geometric scaling as traveling waves 
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We show the relevance of the nonlinear Fisher and Kolmogorov-Petrovsky-Piscounov (KPP) equa- 
tion to the problem of high energy evolution of the QCD amplitudes. We explain how the traveling 
wave solutions of this equation are related to geometric scaling, a phenomenon observed in deep- 
inelastic scattering experiments. Geometric scaling is for the first time shown to result from an 
exact solution of nonlinear QCD evolution equations. Using general results on the KPP equation, 
we compute the velocity of the wave front, which gives the full high energy dependence of the 
saturation scale. 
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1. Geometric scaling p] is an interesting phenomenological feature of high energy deep-inelastic scattering (DIS). It 
is expressed as a scaling property of the virtual photon-proton cross section, namely 

a^(Y,Q)=a^(-^-) , (1) 
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where Q is the virtuality of the photon, Y the total rapidity and Q s an increasing function of Y called the saturation 
scale 0. 

It is convenient to work within the QCD dipole picture of DIS 3J. In the leading logs (Y) approximation of 
Qh' perturbative QCD, the cross section factorizes as 

! r°° r 1 

o 1 (Y,Q)= x Q1 dx 01 dz\rJj{z,x Q1 Q)\ 2 N(Y,x 01 ) . (2) 
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ip(z, xoiQ) is the photon wave function on a qq dipole of size xoi, and z is the longitudinal momentum fraction of the 
photon carried by the quark. N(Y, Xox) is the dipole forward scattering amplitude. Let us define 

^ . 

N{Y,k) = ^ ^lj (kx 01 )N(Y,x 01 ) . (3) 
t— I 1 Jo x oi 

on : 

i In this picture, the geometric scaling property reads 
Oh' 

Within suitable approximations (large N c , summation of fan diagrams, spatial homogeneity) and starting from the 
Oh. Balitsky-Kovchegov equation ,4] , it has been shown ,5] that this quantity obeys the nonlinear evolution equation 

d Y M = a X {~d L )N~aN 2 , (5) 

where a — a s N c /w, x(l) — 20(1) — 0(7) — 0(1— 7) is the characteristic function of the BFKL kernel 0, L — log(fc 2 //c§) 
and fco is some fixed low momentum scale. It is well-known that the BFKL kernel can be expanded to second order 
around 7=5, if one sticks to the kinematical regime 8aY 3> L. We expect this commonly used approximation to 
remain valid for the full nonlinear equation The latter boils down to a parabolic nonlinear partial derivative 
equation: 

d Y M = a X {~d L )N ~aN 2 , (6) 

with 

x(-d L )= X (i) + ^(d L + l) 2 ■ (7) 
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The question we want to address is whether there are exact asymptotic (in Y) solutions of Eq. © exhibiting geometric 
scaling. 



2. Introducing the notation ui — x{\)> D — and defining 7 = 1 — A^/T+8k>/Z), the change of variables 

otD -o , / aZ3 \ 

t=— (l-l) 2 Y , ar=(l- 7 )(L + — y) , 

(8) 



u(t,x) — —— =To 
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D(l- 7 ) 2 Va^li-T) 2 ' 1-7 (!~7) 2 
brings Eq. © for jV into the Fisher and Kolmogorov-Petrovsky-Piscounov (KPP) equation for u: 

d t u(t, x) = d 2 u(t, x) + u(t, x)(l - u(t, x)) . (9) 

We are going to use some general results on the existence of traveling wave solutions to the KPP equation |8j. If 
one chooses an initial condition at time t = to such that u(to, x) decreases smoothly from 1 to as a; goes from — oo 
to +oo, and has the asymptotic behavior 

u(t ,x) ~ ^ , (10) 

x — >+oo 

it is proven Q that the KPP equation admits traveling wave solutions at large times. This means that there exists a 
function of one variable w such that 

V,(t,x) ~ w(x — 771/3 (t)) (11) 
t — >+oo 

uniformly in x. Such a solution is depicted on Fig.Q. The function mp(t) depends on the initial condition: 

mp(t) = c{(3)t + 0(1) for < C , 

m p {t) = 2t- |logt + 0(l) iox0 = c , (12) 
mp{t) = It- flogt + 0(l) iox 0>0 c , 

where c(0) = + 1/0, and the critical value C = 1 corresponds to the minimum of c{0). Note that for critical or 
supercritical > C , the velocity of the traveling wave front dmp(t) / dt is independent of up to a constant. 

3. Let us show that the initial conditions for QCD are of the form of Ea. ljl()(l . For an initial value Y = Yq, corresponding 
to t = to = aD(l — 7) 2 lo/2 and for large L, corresponding to large x, the behaviour of N is given by perturbative 
QCD. Asymptotically in L (i.e. in x at fixed t), it reads 

TV oc k- 2 t" , 70 = 1 . (13) 

We get the relationship between and 7 from the coordinate map ©: = 7o/(l — 7)- In particular, u(to, +00) oc 
Af(Yo,L — » +00) = 0. Formula 113fl through the Fourier transform J3J, corresponds to QCD color transparency, 
namely N(Y, xqi) oc Xqi- 

The condition u(to, — 00) = 1 corresponds to the normalization condition Af(Yo,L — » —00) = D(l— 7) 2 /2 when 
k/ko — > 0. This infrared limit is not attainable in perturbative QCD, but it is reasonable to expect that some upper 
bound on M exist because of strong absorption. 

Interestingly enough, happens to be equal to 2/-v/l + 8u>/D = 1.55276..., i.e. it lies within the supercritical 
regime > C of the KPP equation. Note that in this regime, possible logarithmic prefactors in Eci. (|13fl would not 
modify the traveling wave solution, at variance with the critical case = C ■ Hence color transparency is the only 
relevant property for defining the initial condition at large L. 

According to the above-mentioned mathematical framework, traveling wave solutions to Eq.© for Af exist at 
large Y. We shall show that they correspond to specific geometric scaling solutions. Indeed, the velocity of the 
traveling wave front is given by the third equation in <|12f) because we are in the supercritical regime. This means 
that the large time solutions depend on a single variable x — 2t + % logt. From the coordinate mapping JHJ , one sees 
that we may write the solution as 

M{Y, k)=M f g^py ) wh ere Q 2 S {Y) = k 2 Y~^e aD ^-^ Y . (14) 



u (t, x) 




FIG. 1: Typical traveling wave solution. The function u(t,x) is represented for three different times. The wave front connecting 
the regions u — 1 and u = travels from the left to the right as t increases. That illustrates how the "strong absorption" region 
invades the "transparency" region (see text). 



We have absorbed all 0(1) constants into feg. Ea. ll4|l is nothing else than geometric scaling, Q S (Y) being the 
saturation scale, cf. Eq.|@}. 

4. Equation (|14fl is the exact asymptotic solution of the nonlinear Eq.© and thus geometric scaling appears as a 
universal property of this kind of equations. We obtain the saturation scale including the complete prefactor up to 
factors of order 0(1). It comes from a non-trivial mathematical property of nonlinear equations. To our knowledge, 
it is the first exact result taking into full account the nonlinearities of a QCD evolution equation . 

Comparing our results to previous studies, the exponential growth of the saturation scale with the rapidity Y 
matches the estimates obtained from the extrapolation of the solution of the linear BFKL equation up to the saturation 
scale Q. The powerlike prefactor that we get and prove to be universal,matches the one obtained in Ref.^^j where 
absorption was treated as a boundary condition on the linear equation [llj . 

Several numerical studies of Eq.© have been performed. Evidence for "soliton-like" solutions [l^| and for geometric 
scaling has been found 

[HE1. Recently, formula JT3J was proved to be very well verified |l4j |. This matches the 
mathematical properties we have discussed. 

Among pending mathematical questions, it would be interesting to have a similar approach for the original equa- 
tion J5|, beyond the approximation J7J. The question is whether geometric scaling is still valid in this case or if 
there are scaling violations which could also be universal. Another intriguing question is the meaning of the strong 
absorption bound in the initial condition. We think that the mathematical apparatus for nonlinear partial derivative 
equations could provide interesting hints for these questions. 

As for the physical prospects of our approach, it is inspiring to note that the same equation appears in many 
different contexts, among which the problem of polymers on disordered trees |l5j . In that context, (3 plays the role 
of a temperature and there is a phase transition at /3 = C . (3 > (3 C corresponds to a phase of spin glass type. Similar 
features have been discussed in multiparticle collisions |1(t| . 

R.P. wishes to thank B. Derrida for introducing the subject of traveling waves to him, and E. Iancu for interesting 
discussions on saturation. 
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